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Abstract 



We propose the use of algebras of generalized functions for the analysis of certain highly 
I singular problems in the calculus of variations. After a general study of extremal problems 

. on open subsets of Euclidean space in this setting we introduce the first and second variation 

of a variational problem. We then derive necessary (Euler-Lagrange equations) and suffi- 
cient conditions for extremals. The concept of association is used to obtain connections to a 

■ distributional description of singular variational problems. We study variational symmetries 
and derive an appropriate version of Nother's theorem. Finally, a number of applications to 
geometry, mechanics, elastostatics and elastodynamics are presented. 
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tlT^ ■ 1 Introduction 

o ■ 



The study of singular problems in the calculus of variations has a long history and a sizeable 
literature on diverse aspects of this topic is available (cf., e.g., [T21[nillS] and the literature cited 
therein). In this paper we introduce an approach to variational problems involving singularities 
, based on the theory of algebras of generalized functions in the sense of Colombeau ([H [S]). We 

' are interested in extremizing functionals which are either singular (e.g., distributional) themselves 

or whose set of admissible variations includes generalized functions. 

When considering singular variational problems in the above sense, distribution theory is only of 
limited use due to the nonlinear nature of the typical variational problems. The theory of algebras 
of generalized functions, on the other hand, provides a nonlinear extension of distribution theory 
which allows to model nonlinear singular problems while at the same time providing optimal 
consistency properties with respect to the linear theory. It has found an increasing number of 
applications in linear and nonlinear partial differential equations (cf., e.g., [36] . or [35j for a recent 
survey), regularity theory and microlocal analysis (e.g., [101 [THl [7]), as well as in non-smooth 
differential geometry (e.g., [ini UHl [Ml 131] ) ■ As a rather novel development, the order structure 
in Colombeau type algebras of generalized functions has been investigated in [TJ [301 [32] . It allows 
the formulation of variational problems in the generalized setting and will therefore be crucial to 
this work. 
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Furthermore, we study symmetry properties of variational problems, thereby continuing our in- 
vestigations on Lie symmetries of differential equations in generalized functions ( [Ml [SI [231 HI] ) . 
Our aim here is to derive infinitesimal criteria and establish an appropriate version of Nother's 
theorem. 

The paper is organized as follows. After recalling some standard notations from symmetry group 
analysis and algebras of generalized functions, in Section [2| we examine extremal values of Co- 
lombeau generalized functions on open subsets of Euclidean space. As in the classical calculus of 
variations this will provide the basis for deriving necessary and sufficient conditions for general 
variational problems later on. Already at this stage, new effects occur due to the structural prop- 
erties of the set of scalars in Colombeau theory: they form an ordered ring but not a field. The 
presence of zero divisors requires a more refined analysis than in the smooth case. In Section [31 
we turn to calculus of variations proper. After singling out the appropriate spaces of generalized 
functions for formulating the problem, we study necessary and sufficient conditions for minimiz- 
ers. In particular, we prove a fundamental lemma of the calculus of variations and obtain the 
Euler-Lagrange equations as necessary conditions. We also introduce the concept of minimizer in 
the sense of associations which allows to study variational problems on the distributional level. 
In section [4| we investigate variational symmetries in the framework of algebras of generalized 
functions. Here we derive an infinitesimal criterion and establish a Nother theorem, allowing to 
derive conservation laws from variational symmetries. Finally, section [5] is devoted to a variety of 
applications, focussing on variational problems with singularities that require formulation in the 
Colombeau setting. The first example treats geodesies in a generalized Riemannian metric. Then 
we turn to particle mechanics with singular potentials. The third type of example is concerned with 
elasticity theory, either with degenerate material properties (loss of ellipticity) or with singular 
potentials: membranes with springs attached, beams with discontinuous or vanishing coefficients 
and rods with generalized stress-strain relation. In the last example we turn to the hyperbolic 
case: the wave equation with singular and nonlinear potential. 

We now briefly introduce the notations that will be used throughout this paper. Our basic reference 
for the theory of algebras of generalized functions is [16] . For symmetry group analysis we refer 
to [42]. Concerning the calculus of variations, any standard text on the subject will cover what is 
needed in our approach (e.g., [T^[T31[^ ) 

The spaces of independent and dependent variables will be MP and respectively, or their subsets 
n and U (we set M = QxU). By (x, m^")) we shall denote a point in the n-th order prolongation 
or jet space M^"), whose components are a;*, i — l,...,p, of x and Uj, a = l,...,q, J — 
(ji, . . . , jfe), 1 < jk < P, 1 < k < n, representing u and all partial derivatives of u up to order 
n. The n-th prolongation of a function f : ^ U will be the function pr^"-*/ : — > f/'^"') 
with components representing / and its partial derivatives up to order n. In general, the n-th 
prolongation of any object will be denoted by pr*^"^. By G we shall denote a one-parameter 
Lie group of transformations which acts on M and write for the action (x' ,u') = 5^ • {x^u) = 
^{r],{x,u)) — (2(77, (x, u)), ^'(77(0;, u))). To simphfy notation we will also write <I>^(x,u), S^(x,u) 
and 'i>jj{x,u). G is called projectable if S^(x,u) — S^(x). The infinitesimal generator of G is a 
vector field on M denoted by v and written in the form X^f^i ^'(^' + Sq=i '^al^^''")^""- 
The total derivative with respect to will be denoted by Di. 

The algebras of generalized functions on which our approach is based will be those of special (or 
simplified) type ([13 Ch. 1]), so we will henceforth drop the adjective "special" and simply refer 
to them as Colombeau algebras. The algebra G{^) is the quotient Sm{^)I-N'{VI)^ where 

£M{n) = {(ue)e e (C°°(17))(°'il|Vi\: CC O Va e N[; e N : sup \d°'u,{x)\ = ©(e"^) as e ^ 0} 
Af{n) = {(ue)e e C°^{nY°'^^\yK CCnVae N[;Vm e N : sup \d°'ue{x)\ = 0(£'") as £ ^ 0}. 

J\f{n) is an ideal in the algebra £]^J(fl) (all operations are defined componentwise, i.e., for fixed 
e). Elements of £m{^) and J\f{V,) are called moderate, resp. neghgible nets of smooth functions. 
C/(r2) is an associative, commutative, differential algebra whose elements are equivalence classes 
u = [{ue)e]- ^ ^ G{^) is a fine sheaf of differential algebras on Rp. In particular, there is a 
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well-defined notion of support in G(Cl) and we denote by Gd^) the space of compactly supported 
generalized functions in Q{il). Moreover, G{^) contains the space as a linear subspace and 

C°°{fl) is a faithful subalgebra: embedding is effected basically by convolution with a fixed moUifier. 
This embedding is a sheaf morphism that commutes with partial differentiation. The concept 
of association provides a means of assigning macroscopic properties to elements of Colombeau 
algebras: elements u, v of a Colombeau algebra are called associated if — ^ in 2?'. u is 
associated to the distribution w if — s- w weakly. These notions are independent of the chosen 
representatives. 

Given another open set il' C M?, we may consider those elements of Q{fl)'^ possessing a represen- 
tative such that Ue(ri) C O' for all e and which is compactly bounded or c-boundcd in the 
sense that 

VA' CC n 3K' CC n' 3ea e (0, 1] Ve < Eq : MK) C K' . 

The space of all c-bounded generalized functions from fl to il' is denoted by Elements 
of Q[n,il'] can be composed unrestrictedly as well as inserted into elements of G{^')- A similar 
definition can be given for smooth manifolds instead of open sets and a functorial theory of 
generalized functions based on this notion has been developed in [25l |30] . 

The algebra of tempered generalized functions Gt{^) is defined as the quotient £r{i^)/AfT{^), 
where 

£^{n) = {{ue)e e C°°(r!)("^il|V« e 37V e N : sup(l + \x\)-^\d"u,{x)\ = 0(£"^) as e ^ 0} 
Afr{n) ^ {{u,), e C°°(f7)("^il|Vae N[; 3/ e N V™ G N : sup{l + \x\)-^\u,{x)\ = 0{e"') as e ^ 0} 

The algebra Gt can be used to implement Fourier transformation in the Colombeau setting (anal- 
ogous to the space S' in distribution theory, which in turn is embedded into Gt)- Our main 
interest in it, however, relies on the fact that insertion of elements of G into elements of Gt yields 
well-defined elements of G- 

We shall also make use of the following 'mixed' variant of Colombeau algebra: let Gt{^ x = 
Srifl X n')/Afr{n X rj'), where n' is a subset of M^', and 

£r{n X n') = {{ue)e G <c°° {n X n'Y°'^^\yK ccnyae 3n en-. 

sup (1 + \y\)-''\\d^u,ix,y)\ = 0(£-^) as £ ^ 0} 

xGK.y£fl' 

Afri^ X n') = {{ue)e G {il X n'Y"'^^\yK CC f7 Va e Ng+P' 31 Vm e N : 

sup {l + \y\)-'\\d"u,{x,y)\ = 0{e"') as e ^ 0} . 
xeK.yen' 

Thus, its elements satisfy tj-estimates in the il-variable and ^T-estimates in the fi'-variable. 
It may be seen as a nonstandard feature of Colombeau algebras that, due to the presence of 
infinitesimals, their elements are not uniquely determined by their values on the usual points in 
their respective domains. In order to achieve such a unique determination one has to use so-called 
generalized points in the following sense ([22]): The space of generalized points f2 is the set of all 
equivalence classes in 

= {{xe)e e f^^"^^l|3iV e N : \xe\ = 0{e-^) as e 0}, 

where the equivalence relation is defined by 

(xe) - {ye)e ^ Vm G N : - y,\ = 0(e") as e 0, 

i.e. f2 = fli// ~- The set of compactly supported points is 

[{x,),] e Cl\3K CC n 3eo e (0, 1] Ve < eo : a;e e K}. 

If = RP we shall write W resp. M.^. As a special case, when $7 = M, we obtain E, the ring 
of constants in any of the above algebras. Elements of G{^) are uniquely determined by their 
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point values in all compactly supported generalized points in fi. If is an n-dimensional box and 
u € Gt{^) then 7i = in Gt{^) if and only if u{x) = in R, for all x G (l. Moreover, one can 
show that Colombeau generalized functions are in fact uniquely determined by their values in all 
near-standard points. Here, a point i e is called near-standard if there exists x G Q such that 
x^ ^ X as e ^ for any representative (xg)^ of x. 

The above notions can be employed to define generalized Lie group actions in the Colombeau 
setting (cf. [Ml [3 ES] ) : A generalized group action on W is an element $ e '?r(]R x IR^)^ which is 
a one-parameter group in the sense that $(0, . ) = id in C/^(RP)p and $ (771 -I- 772, •) = ^{rji, ^{112, ■ )) 
in t/r(M^ X MP)P. If ^ G QrOSP)P is a generalized vector field with the property that there is a 
unique generalized group action $ e GtO^ x W)p satisfying ^$(77, x) = ^(<i>(?7, a;)) in ^^(R^+J')^' 
(i.e., if <f> is the unique flow of ^) then ^ is called the infinitesimal generator of $ and both ^ and 
$ are called tj-complete. If, in addition, <I> and ^ have representatives {^e)e and (^e)e such that 

is the flow of for each e G (0, 1], then <I> and ^ are called strictly tj-complete. 
The remainder of this introduction will be devoted to the notion of positivity in the generalized 
functions setting. Since the results in this area are rather recent and come from a variety of 
different sources, we collect below those parts of the theory which will be crucial for our further 
investigation (cf . [H SOI [32] ) . 

By [16, Theorem 1.2.38], an element x of M is invertible if and only if it is strictly nonzero in the 
following sense: for any representative (xe)^ of x there exist a > and sq such that \xs\ > e", for 
all £ < eg. This, in turn, is equivalent to x not being a zero divisor Theorem 1.2.39]). 
i € R is said to be nonnegative, 5: > 0, if there exists a representative such that each Xe is 
nonnegative. Equivalently, for each representative (a;e)e of x and each a > there exists Eq such 
that Xe + s°' > 0, for all e < eg- i G M is said to be strictly positive, 5; > 0, if it is positive and 
invertible. By the above, this means that there exists a > such that for each representative (xe)^ 
of X there exists Eq with x^ > e°, for all e < Eq. Note that since R contains zero divisors, a: > is 
strictly stronger than i > and i ^ 0. We shall need the following characterizations from [3^: 
i e R is strictly positive if and only if 

(i) X is invertible and has a representative (x^)^ such that x^ > 0,Ve G (0, 1]; 

(ii) for each representative (a;e)e of x there exists an Eq such that x^ > 0,Ve < £9. 
As a corollary we have: if a; is a near-standard point associated with a; G R then 

(i) if X 7^ (resp. a; > 0) then x is invertible (resp. strictly positive); 

(ii) if i > then a; > 0. 

Note that the converse implications do not hold as demonstrated by the following simple examples: 
take X :— [(£)e]- Then x is invertible (moreover strictly positive), but a: 0. Also, [(— £)e] « 0, 
but [(— £)e] < 0. We shall repeatedly make use of the following result: 

Lemma 1.1 If x,y G R and \x\ < s ■ \y\, for all invertible s G R with < s < sq for some 
< So G R, then i = m R. 

Proof. Suppose a; 7^ in R. Then for any representative (a:^)^ of x there is Z G R and a sequence 
£fc \ such that 

\x,J>Ei, V/cgN. 

Further, by assumption, for any representatives (jje)^ of y and (sg)^ of s and for all a > 

3£o V£ < £0 : s,- \y,\ -\x,\+e''> 0. 

In particular, this inequality holds for each Ek and /c G N large enough. Fix {ye)e and choose 
AT G N and /cq G N for which \y^^ \ < £fe ^, V/c > fco. Then 

S • > £' - £" > f'+I 
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for a, k large enough, hence 

> ^N+l+l 

Since this inequahty would have to hold for all invertible s G M, < s < sq, we arrive at a 
contradiction (e.g. take s = [{£'^~^''~^^)e])- ^ 



Remark 1.2 The above proof in fact shows that it suffices to check the assumption for all s = 
[{e"')e] with m > 0. 

Turning now to non-constant generalized functions we recall from [IDl the definitions of non- 
negative and strictly positive generalized functions as well as a characterization of positivity: 
/ e g{n) {n C W open) is said to be nonnegative, / > 0, if V {fe)e CCl7 Va > Beq Ve < 
Eq : ini^^K fs{x) -l- > 0. / is said to be strictly positive, / > 0, if V(/j)j CC Q 3a > 
3eo Ve < £o : inix^K fe{x) > e°. The latter can be characterized as follows: a generalized func- 
tion / e Q{Q) is strictly positive if and only if V(/e)e CC fl 3eo Va; G Ve < eo : /e(x) > 0. 
Moreover, by [16} Theorem 1.2.5], / is invertible in G{^) if and only if |/| is strictly positive. 
We conclude this brief overview by recalling the notion of positive (semi)definiteness of symmetric 
bilinear forms on Rp. For this purpose we first need the definition of a free generalized vector (cf. 
[32]): let I G MP and I ^ in RP. We call f free if whenever A G K and A • f = then it follows that 
A = in M. Let A G Rp be a symmetric matrix (i.e. A = in Rp ) . A is positive semidefinite if 
for any ^ G Rp, C^A^ > in R. A is positive definite if for all free vectors ^ G Rp, ^^A^ is strictly 
positive (i.e., invertible and positive) in R. Equivalent characterizations for positive definiteness, 
resp. semidefiniteness of a symmetric matrix A G Rp are: 

(i) all eigenvalues of A are strictly positive (resp. nonnegative); 

(ii) A has a representative (^e)e such that each is symmetric and positive definite (resp. 
positive semidefinite); 

(iii) for each symmetric representative {A,,)^ of A there exists Eq such that Ag, is positive definite 
(resp. positive semidefinite), for all e < Sq. 



2 Differential calculus in algebras of generalized functions 

Basic necessary and sufficient conditions for extremals in the calculus of variations arc usually 
based on the following elementary results on real valued functions: Let / map Q to R, where is 
an open subset of RP. Let a;o G be a local minimum (or maximum) of / G C^(Sl). Then all first 
order partial derivatives of / vanish at xq. Such points are called critical (or stationary) points. 
When / G C^(r2) then, in addition, the Hessian matrix of / is positive (or negative) semidefinite 
a.t Xq. a sufficient condition for a local extremum of / G C^{il) is as follows: if xq G is a critical 
point of / such that D^f{xo) is positive (or negative) definite on Rp, then xg is a local minimum 
(or maximum) of /. 

In this section we want to derive analogous criteria for extremals of generalized functions. The 
main difficulty in carrying out this task lies in the algebraic structure of R: the existence of 
zero divisors necessitates certain adaptations in the classical arguments and provides for a more 
complex structure of the attainable results. In this paper we are interested mainly in the extrema 
of generalized functions in classical points (i.e. in a; G Rp) as this suffices for the applications in 
the calculus of variations we aim at. We will however add remarks which give analogous criteria 
for extrema in generalized points. 

As in the classical case, we define a local minimum (resp. maximum) of / G G{Q) (fi is an open 
subset of RP) to be an G with the property that there exists a neighborhood il' C of xq 
such that f{xo) < f{x) (resp. f{x) < /(xq)), for all x G i^^- 
Our first result provides necessary conditions for a minimum: 

Proposition 2.1 Let f G G(fl). If xq E is a local minimum of f then 
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(i) Df{xo) = m W; 

(a) D^f{xo) £ MP is positive semidefinite. 
Proof. Without loss of generality we may suppose that xq ^ 0. 

(i) For any i £ {l,...,p}, let gi{-) :— /(O, . . . , 0, •, 0, . . . , 0) be the i-th partial function. Then 
gi £ G{^i) i^i = {t £ R \ {0, ... ,t, ... ,0) £ ft} open in M) and is a local minimum of gi on 
Hi. We have to show that g[{Q) = 9^/(0) = in R for all i. Let i £ {1, . . . ,p} and assume that 
minimizes gi over (— so,so)c (so > 0). Let s = [{se)e] £ (— so,so)c ^^.^ take any representative 
{gie)e of gi. Then for all e we have a Taylor expansion of gi^ at 0: 

with < 9^ < 1. gi{0) < gi{s) says that 

Va > 3eo Ve < £o : fc(se) - g^eiO) + > 0. 
Set Me := sup,g[_,„,,^] \gUx)l Then 

< Se5L(0) + 'fgUOeSe) + s'' < s.gUO) + ^Af, + e^ (1) 

Suppose s is strictly positive (hence invertiblc). Then for any representative {se)e of * there exist 
m £ N and ei such that > e™, for all e < ei. Divide ([T]) by for e < min(£o, 

0<.4(0) + ^ke|+e^ 

where a' = a — m. Otherwise, if [(se)e] is strictly negative we similarly obtain 

M , 
Q>g[M~^\-^e\-e'^ ■ 

Collecting both cases yields 

for all invertible generalized numbers s = [{s^);,] £ (~sq, sq)^, with M ~ [(M^)^]. The claim now 
follows from Lemma [Ol 

(ii) Taylor expansion of / around gives: 

Mx,) = A(0) + DfMM + ixj • D'MO) ■ X, + ^iix,Dff,mx,). 

Now for any a > 0, fe{xe) - /e(0) + e° > for £ sufficiently smah. Also by (i), Df^{0){x^) = 
0{e°'), for all a > 0. Moderateness of D^f yields the existence of some > such that 
suP|3^|<i \\D"^f{y)\\ — 0{e~^). Thus if we suppose that \xe\ = with s = ^^i^ then we ob- 
tain that xj ■ D^f^{0) ■ Xe + e""^ > for any given a if e is sufficiently small. 
D^f{0) possesses a symmetric representative, so for each e there exists an orthogonal matrix 
such that with U [iUe)e] £ U-D'^ f{0}U'^ = diag(Ai, . . . , Ap) with Ai, . . . , Ap the eigenvalues 
of 1)^/(0) in K (cf. [32l Lemma 1.6]). Denote by x\. the corresponding unit eigenvectors (the 
columns of U) and set yi^ := e'^ ■ x\^^. Then by the above, for each a' > and e small 

< + yl ■ D^fM ■ y,, = e'^' + e^^A.^ , 

implying that all A^ are positive in M. As mentioned in the introduction, this is equivalent to 
positive semidcfiniteness of D^f{0). □ 
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Remark 2.2 {Extrema in generalized points). In order to extend the validity of 12.11 to aio G i^c 
we first need a concept of local minimum in this setting. We employ the notion of support of a 
generalized point xq G fic introduced in [9 j the support supp(5;o) of io is defined as the set of 
accumulation points of any representative (a;oe)e. We call Xq a local minimum of / G Q{^) if there 
exists a neighborhood fi' of supp(a::o) such that /(io) < f{x) for all x € Cl[. The proof of 12.11 (i) 
can be adapted to this situation: set gie{t) '■= fe{x\^, . . . , x^^ + t, . . . , x^^). Then although there 
need not be a finite interval on which all gi^ are simultaneously defined, the above arguments may 
still be applied for e small fusing [L2|l and each i, showing that D/^xq) = 0. Similarly we may 
modify the proof of 12.11 (ii). 

The above result demonstrates that the typical necessary conditions for local minima are entirely 
analogous to those for minima of classical functions: the gradient vanishes and the Hessian matrix 
is positive semidefinite. However, as the following example shows one should not expect such 
a direct analogy with the smooth setting in the case of sufficient conditions: We construct a 
generalized function whose first derivative vanishes at zero and whose second derivative at zero is 
strictly positive but which does not have a minimum in 0. 

Example 2.3 Let ip e C(f (M), supp^s C [-2,2], (^|(_i/2.i/2) = x"^ and (^(-1) = (^(1) = -1. 
Define /,(•) := ^{-) and / := [(/,),]. Clearly, / G g{R), ^(0) = 0, /^(O) = and /^'(O) = ^, so 
/"(O) is strictly positive (even /"(O) oo). But zero is not a local minimum since the value of / 
in the generahzed point [(e)e] is —1. 

We therefore have to alter the condition of positive (semi)definiteness in zero to also take into 
account neighboring generalized points. Of course we would like to include as few points as 
possible. The following example shows that merely considering near-standard points associated 
to is insufficient: 

Example 2.4 Let ip G C^(R), supp(^ C (—1, 1) and </3|[-i/2.i/2] = ^1- E)efine 

oo , 
n=l ^ e ^ 

and F := [{F,)^]. Then F G ^(M), F(0) = and F'(0) = in M. Moreover, F"{i) > in R, for aU 
near-standard points 5; « (in fact F"{x) = in R for all x Ki 0). But F does not attain a local 
minimum at zero. For any no G N, F^{-^) ~ — e"" for 2e < I/uq — l/(no -|- 1), hence F{:^) < 
in R. Therefore, is not a local minimum of /. 

As suggested by the previous examples we therefore give a sufficient condition for a local minimum 
which is based on the behavior of / in a neighborhood of the critical point. 

Proposition 2.5 Let f G ^(^1), xq G il, Df{xo) = and let D^f{x) be positive semidefinite for 
all i G (with f2' a star-shaped neighborhood of xq in Q). Then xq is a minimum of f in Q'^. If 
in addition D'^ f(x) is positive definite for all i G then the minimum is unique on Q,'^. 

Proof. Again, taking = is no restriction. Let x = [{x^)^] G Vl'^. Take a representative {fe)e 
of /. For each e we have 

fe{x,) = A(o) + DfM{^e) + \{{x,Dff,){e,x,) 

with G (0, 1). Since D/(0) = in R it follows that for aU m > and e small, \Df,{Q){xe)\ < 
e™. On the other hand, by the assumption that D^f\(y is positive semidefinite we know that 

Va' > 3eo such that l/2{{xi,Dy f^){9^Xi.) + > 0, Ve < eo- Let a > be given and set a' and 
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TO to be a + 1 . Then 



< e''+' + Mx,)~MO)-DfM{xe) 

< e-+' + Uxe)-feiO) + \DfMiXe)\ 

< £'^+i + /e(a;e)-A(0)+e'^+i 

when £ ^ 0. Hence f{x) — /(O) > and the first part of lemma is proved. 

To show uniqueness, let x ^ Cl'^ and i ^ in MP. Then there exist I and Sk \ such that 
l^Ekl > ^fc' for fc e N. Choose a representative {fe)e of / with Df^{0) — for all e and write 
fe{xe) — /^(O) + ^{{x^D)^ f^){9sX^). In order to bring to bear the positive definiteness of D^f we 
embed the sequence {xei^)k into a free vector in il^ as follows: Let y be any free vector in fi^ (e.g., 
y = [(e)e]) and define x' as [(a;g)£], where 




y^ , otherwise 



Then x' is free, so the positive definiteness of D^f in z :— [{6^Xe)e] G ^'c implies strict positivity 
of {x')~^ ■ D'^ f{z) ■ x' . Thus there exists some m > such that for large enough fc g N 

This implies that f{x) ^ /(O) in R, for all x £ il'^, and uniqueness follows. □ 



Remark 2.6 As in 12.21 the above result may be extended to generalized points xo, if we suppose 
il' to be a convex open neighborhood of supp(io)- 

The uniqueness condition in the previous result raises the question whether positive definiteness 
of the second derivative could be replaced by the assumption of D^f being positive semidefinite 
but different from zero in M. These two notions are different in R, due to the existence of zero 
divisors. The following example, however, shows that the latter condition is too weak. 

Example 2.7 For any zero divisor a > in R the assignment 

f{x) := ax'^, 

defines an element of Q(R) which attains its (even global) minimum in the point 0. Moreover, 
/'(O) = in R and /" = 2a is positive and different from zero everywhere in R. Nevertheless the 
minimum at is not imique: take any s 7^ in Rc with a ■ s = 0, then also /(s) = 0. 

We conclude this section with a characterization of nonnegativity, resp. strict positivity, of gen- 
eralized functions by their values in all near-standard points. As a consequence we note that the 
respective assumptions in the previous results can equivalently formulated in terms of this proper 
subset of generalized points. 

Lemma 2.8 Let f G G{^)- Then / > m G{^) (resp. f > in G{^)) if and only if f{x) > in 
R (resp. f{x) > inR) for all near- standard points x G flc- 

Proof. We give a proof for the characterization of nonnegativity of /, the case of strict positivity 
can be established along the same lines. The condition is clearly necessary. Conversely, suppose 
that / is not nonnegative in G{fl). Then 

3(/e)e 3a>0 3K CCnykeN 3ek < j 3xk e K : feA^k) + < 0. 
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Without loss of generality we may assume that Sk+i < £k and that the sequence {xk}keti converges 
to some xq G K. Set :— Xk for e g (£^+1,6^] and x :— [{xe)e]- Then x is a near-standard point 
and f^i^{xe^) < Vfc e N, hence f{x) ^ 0, a contradiction. □ 

Inspection of the above proof shows that in fact the subset of so-called sequential near-standard 
points suffices to characterize nonnegativity resp. positivity of generalized functions. Here, a near- 
standard point X ~ X is called sequential if it has a representative (xe)^ of the following form: 

3s k \ 3Xei^,Xek X, SUCh that Xe = Xei^ye e (£fc+i,efc]. 



3 First and second variation 

Let n be an open subset of R" and let D ~ fl or D — Cl. Let be a subset of G{D) and let 
£ : ^ — > M. We want to determine extremals of C under certain admissible variations. To make 
this problem accessible to an analytical treatment we make the following assumptions (which will 
automatically be satisfied in the setting most relevant to us, namely ([3]) below). Let the space 
of admissible variations be an M-submodule of G{D) such that A + Aq C A. Moreover, we suppose 
that for all It € and all v G the mapping 

is an element of G{J^'^), where J"'" is some open interval around in M. 

We call u & A a minimizer of the functional £ in ^ if C{u) < C{w) in K for all w G A. Finally, 
we define the first and second variation of £ at u G .A in the direction v G Ao as 

C-^^s) (2) 

Under these assumptions we have 

Proposition 3.1 Let u A be a minimizer of the functional C in A. Then 

(i) 6C{u; v) — inM. for all v G Aq; 

(n) S'^C{u; v)>0 inR for all v G Aq- 
Proof. This is immediate from Proposition [Til ^ 

Turning now to our main object of study, let Q be an open, connected and bounded subset of W 
with smooth boundary dft. Let £ : t/(ri)^ — > R be a functional of the form 

C{u) := [ L{x,u'-''Hx)) dx, (3) 
Jo. 

where L G Qt{{^ x R'^)'-"-'). We will write for the functional u ^ f^-^L^{x,u) dx, where 
u G C°°{Cl) and {L^)^ is any representative of L. 

We seek a minimizer u among all functions in G{fiy which additionally satisfy the boundary 
condition = "o, uq G Gidfl)"^ (here dfl is viewed as a (p— l)-dimensional submanifold of Rp). 
The integrand L is called the Lagrangian of the variational problem. In the notation introduced 
above we have 

A = {ivegin)'^\w\9n=uo} (4) 
for some uq G Q{d^)'^, with the corresponding space of admissible variations given by 

Ao = {weginy\w\an = 0}. (5) 

Note that the integral in ^ is well-defined since the composition of a tempered generalized function 
with an element of C/ is a well-defined element of G (cf. [ini Proposition 1.2.30]). Also, ^ is the 
most general form of the variational problem, involving p independent, q dependent variables and 
all partial derivatives of dependent variables up to order n. In some cases, in order to simplify 
notation we are going to consider only the case of one unknown function; the generalization to 
two or more dimensions is then a technical matter. 



6Ciu;v) - 
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Remark 3.2 In the formulation of the variational problem, instead of tempered generalized func- 
tions one may alternatively use c-bounded generalized functions: assume that L € ^;((f7xM9)(")). 
Since c-bounded generalized functions may be composed unrestrictedly (see [161 Proposition 
1.2.8]), a minimizer then can be sought among all c-bounded generalized functions with pre- 
scribed boundary conditions. In fact, for the integrand to be well-defined not only u but also 
all partial derivatives of u up to order n have to be c-bounded, i.e. the admissible set is A 
{w\ w'"-* e (R"^)*^"'] A w\dn = uo}. Although convenient for an invariant geometric formula- 
tion of the theory, this latter requirement necessitates certain unwanted restrictions on the degree 
of singularity of u (e.g., (5-like singularities are excluded). We will therefore mainly operate in the 
tempered setting, but will simultaneously keep track of the c-bounded case. 

In order to derive Euler-Lagrange equations in our setting, as in the classical case we first derive 
a fundamental lemma of the calculus of variations: 

Lemma 3.3 (Fundamental lemma of the calculus of variations) Let u £ If 

u{x)ip{x) dx = m R 



for all If G Gci^) then u = in G{Q). 

Proof. Suppose to the contrary that u = [{ue)s] with (u^)^ ^ A/'(il). Then 3K CC ft 31 > QVk e 
N 3ek \ 3xk e K such that \ue^^{xk)\ > e^. On the other hand, (ue)e is moderate, so 3N > 
such that sup^g^ \Du^{x)\ < £~^. We may choose some ko such that for k > ko and any x £ Q 
with — Xfel < ^s''/^^ , the entire line segment connecting x and Xk is contained in Q. Then 

Uejx) = Uejxfe) -I- Du^^{xk + Ok{x - Xk))ix - Xk), 

for some 0k £ (0, 1), so that [^^^^(x)! > e^. — £k^\x — Xk\ > ^^k for all x with \x — Xk\ < l^fe^^- 
Choose If £ C^{Bi{Oj),ip > and = 1 on Bi{0) {Br{0) denoting the open ball of radius r 
around 0). Define ip^ £ 'D{il) by 

(pe{x) := e~^~^ Lp{e^^-^ {x - Xe)) ■ Slgn{Ue(Xe)) 

where x^ = x^^ for all e £ (sfe+i, £k\ and all k > k^. Then {p := [(<y5e)e] £ Qc{^)- Since m^j. does not 
change sign on the support of (/s^^ and since u^^{x) ■ Pcki^) ^ h^k^ whenever |x — x^j < 
we obtain 



-B 1 L + N{^k) 



> / u^k{y)Vek{y) dy 



> const. 

which implies that {j^ Ue(y)(pe{y) dy)^ ^ A/'(ri), a contradiction. Therefore, u = in Q{Q). □ 

As in the smooth setting, this result allows to derive the Euler Lagrange equations as necessary 
conditions for a minimizer of ([3]): 

Theorem 3.4 Let u £ A be a minimizer of Then u is a solution of the system of Euler- 

Lagrange equations 

E^{L) = Q, a = l,...,g, ing{^l), (6) 
where Ea is the a-th Euler operator 

(the sum extends over all multi-indices J — (ji, . . . ,jk), ^ < jk < P,^ < k < n). 
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Proof. By Proposition [3TT] (i), 5C{u;v) = 0,Vu G ^o- Now 

SCiu;v) = / 4l L{x,{u + svY''\x)) dx 
Jn asL=o 

Integration by parts gives (with D denoting the total derivative) 



SC{u; v) 



I il[T.^'D)j^{xM-H^))\-v%x)dx 
Jn L ^ auj J 



J 

The claim therefore follows from Lemma [5751 □ 

Remark 3.5 In the case n = 1, i.e. when C{u) ~ Jj^ pr^'^^M(x)) dx, the a-th Euler operator 
takes the well known form: 

i—i '■ 

where uf — dvP^ jdxi. 

The following result provides sufficient conditions for niinimizers of the variational problem ([3]) . 
Proposition 3.6 With L as in \^ and A, Aq as in suppose that 

(i) 5C{u\ w) = m M for all v G Aq. 

(ii) S'^C{w; v) > for all w CI A and v G Aq. 

Then u is a miniraizer for C If, in addition, s i— > 6^C{u + sv; v) ^ in Q{[0, 1]), Vu G A^, v ^ Q, 
then the minimizer is unique. 

Proof. Given any w G A, set 'i' s i-^ C{u + s{w — u)) G By assumption we have 

*'(0) = 5C{u; w - m) = (since w - u Aq) and 4'"(s) = 5'^C{u + S(u; -u);w-u)> 0, VS G Mc- 
We now apply Lemma [2.51 to conclude that is a minimum of 'J, i.e. ^'(0) < 5'(s),Vs G Mc. In 
particular, for s = 1 we have that C{u) < C{w), and since w G A was arbitrary the claim follows. 
If, in addition, s ^ 5'^L{u + s{w - m); w - u) 7^ in Q([0, 1]) for aU w & A, i.e. s 1-^ *"(s) 7^ 
in t/([0, 1]), then according to [T71 Lemma 3.1] (which in fact is closely related to our Lemma [575)1 
we obtain that ^ 

«'(!)- ^'(0)= / {\ ~ s)^" [s) ds ^ Q . 
Jq 

Therefore C{w) ~ C{u) ^ for all w ^ u, and uniqueness follows. □ 

Remark 3.7 The uniqueness criterion in Proposition l3.6l can be adapted to provide an alternative 
to the uniqueness condition in Proposition l2.5[ to wit: under the assumptions given there, suppose 
additionally that s ^ D'^f{xo + s{y — xo))[y — xo,y — xq] 7^ in Q[0, 1] for all y G Then the 
minimum is unique in Q'^. 

Example 3.8 (Quadratic forms). Suppose that for each e G (0,1], ae : C°°{n) x C°°{n) M 
is symmetric, bilinear and continuous. We additionally suppose that {a^{u^,v^))^ G £m for all 
{us)e,{'"e)e ^ ^^/(f^), and {a^{u^,Vs))s G TV for (u^)^ G £m{^) and {v^)^ G 7V(n) (where R = 
£m/J^), so that 

a : g{n) X g{n) m 

a{u,v) := [{ae{Ue,Ve))s] 
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is a well-defined bilinear map. Similarly, for each e E (0, 1] let : C°°{ft) ^ R be linear and 
continuous and suppose that {fg{ue))e & £m (resp. Af) for each (u^)^ G £m{^) (resp. J\f{Ct)). 
Then set / : g{n) ^ t, f{u) := [(/, o u,),]. 
We consider the functional 

£{u) := ^a{u,u) + f{u) . 

Again we will write for the functional u i— > ^a^{u,u) + /e(u). In this situation, condition (i) 
from Proposition 13.11 reads 

SC{u,v) = aiu,v) + fiv) = (8) 

and d'^C{u,v) = a{v,v). Thus if a is positive semidefinite in the sense that a(u,u) > for all 
u e G{(1) then any solution of ([5]) is a minimizer of C. 
As a concrete example let 

C{u)^l- / {a{x)\Wu{x)\^ + I3{x)u{x)^) dx + / -f{x)u{x)dx (9) 

with A, Ao as in (g]), ([5]), a,/3,7 € and a > 0, /? > on fj. Under these assumptions, a 

is positive semidefinite and by the above, any solution of ([8]) is a minimizer of £. Moreover, the 
minimizer is unique by Prop. [3^ in case a is invertible on Q. In fact, if 

S^C{u;v) ^ a{v,v) ^ [ {a{x)\\7v{x)\^ + (}{x)\v{x)\^) dx ^ 



it follows that Vw = in G{^), so w is a generalized constant. If v is an element of Ao, this 
constant must be zero. Therefore, S'^C{u,v) ^ for all v G Aq^ f 7^ 0; this in turn implies the 
hypotheses of Prop. 13.61 

On the other hand, uniqueness may be lost if a is not invertible: let — (—1, 1), and 

£(u) = - y q;|u'(x)P da; + ^ au{x) dx 

with a > a zero divisor, say acu — 0. Then 5'^C{u, v) = J^^ a\v'{x)\'^ dx > for all v. However, 
both u{x) = x'^/2 — 1/2 and i2(a;) = ujx'^ /2 — uj/2 are minimizers, because they satisfy the Euler- 
Lagrange equation 

a( - u"(a;) + 1) = 0, m(-1) = m(1) = 0. 

We now turn to an analysis of variational problems in the distributional setting. In the language 
of Colombeau algebras, this amounts to expressing extremal properties in terms of the concept of 
association: 

Definition 3.9 An element u of G{fl) is called a minimizer in the sense of association for the 
functional ^ if for any (hence all) representatives {Cs)e of C and {ue)e of u we have 

\im{£,{u^ + ip)-£,{ue))>0 VipeV{n). (10) 

£—♦0 

For a large class of quadratic forms, the following result provides a necessary and sufficient condi- 
tion for a minimizer in the sense of association: 

Lemma 3.10 Let C be a functional of the form ^ where a is positive semidefinite and such 
that A{(p,ip) :— lim^^o exists for all (p G 'D{fl). Let u G G{^). Then the following are 

equivalent: 

(i) u is a minimizer in the sense of association for C 
(a) lime^o(ae(w£, (ys) + /e((^)) = for all ip G T){n). 
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Proof. (ii)=>(i): Let ip G I?(fi). Then Ce{u^ + ^) - ^e{ue) — a^{us,ip) + ^ae{(p,ip) + feif) 
iA(^,</.)>0 (e^O). 

(i)=>(ii): Again fix (^9 G ^{Q). Inserting np instead of ip in pop . we obtain 
\im{a,iu„ p) + f,{p)) + l-A{p, p)>Q Vr > , 

e— >0 Z 

and 

lini(a,(w„ ip) + f,{p)) + (^) < Vr < , 

£— »0 Z 

wfiich together give the claim. □ 
The following is a classical example due to Weierstrass, cf. also [TS]: 
Example 3.11 Let 57 = (—1, 1), c 7^ d G K and consider the functional 

C{u) = j x^u'{x)^dx. 

We want to minimize C subject to the boundary conditions u{—l) = c, u(l) = d. Here, is in 
fact independent of e. Denote by v the step function equal to ci on (—1,0) and equal to C2 on 
(0, 1) {v is the "expected" minimizer). Let p be a standard moUificr (p G J p = 1) and set 

Us := V * pe- Then 

\im {C{u, + p) - C{u,)) = lim [ {2x'^u'{x)p'{x)+x'^p'{xf)dx>0 ^p £V{n) , 

so u ~ [{ue)e] is a minimizer in the sense of association. Moreover, C{ui,) for e — )■ and u 
is associated to v. By Lemma [3.101 (or by direct calculation), u solves the corresponding Euler- 
Lagrange equation in the sense of association. 

Note, however, that u is not a minimizer in G{—1, 1). Otherwise by Theorem l3.4l it would have to 
satisfy the corresponding Euler-Lagrange equations (with equality in Q). This, however, only has 
the classical solution, so also in G{—1, 1), no solution satisfying the boundary conditions exists. 

For more general variational problems of the form ^ we need additional conditions in order to 
have Euler-Lagrange equations (as in Lemma l3.10l fii)) as necessary conditions for minimizers (in 
the sense of association): 

Proposition 3.12 For the variational problem suppose that the total second order derivative 
of L w.r.t. the variables u^"' is globally bounded (independently of e). If u £ G{^) is a minimizer 
of C in the sense of association, then u satisfies the Euler-Lagrange equations for C in the sense 
of association: 

Ea{L)^Q, a = 1, ... ,9, 



I.e., 



lim / Ea{L,){ue){x)p{x)dx Vv? G X>(f)) Va . 



Proof. To simplify notations we give the proof for p = q = 7i ~ 1, the general case being entirely 
analogous. By Taylor expansion we have 

Le{x, y + h, z + k) ^ Ls{x, y, z) + d2Ls{x, y, z)h + d^Leix, y, z)k + Re{x, y, z, h, k) 

where ^ 

Reix,y,z,h,k) = ( / dI ^Ls{x,y + th, z + tk){l - t) dt)[{h, k), {h, k)] 
Jo 
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and by assumption, \R^{x,y,z,h,k)\ < C\\{h,k)\\'^ for all e. Let ip G r G R. Then 

< lim(£j(ii£ + T(p) — Ce{u^)) 

= lim ( r / d2Ls{x,Ue{x),u' {x))Lp{x)+d:>,Le{x,Ue{x),u' {x))ip'{x)dx 

+ i Re{x, Ue{x), u'^{x), TLp{x) ,Tip' {x)))[{ip{x) , if' {x)), {ip{x), <^'(x))] dx 

Jn 

Since the second integral is uniformly bounded, dividing by r and distinguishing the cases r > 
and T < (cf. the proof of Lcmma FS-lOp implies the result. □ 



4 Variational symmetries and Nother's Theorem 

In this section we turn our attention to symmetries of variational problems involving generalized 
functions. After a quick recapitulation of the relevant notations and definitions in the smooth 
setting (following [55', Ch. 4]), our goal is to introduce an appropriate notion of generalized varia- 
tional symmetry groups, retaining the fundamental properties of classical variational symmetries. 
Based on this, we will then derive a corresponding version of Nother's theorem. 
Turning first to the smooth setting, let C : C°°{il) ^ R be a functional of the form 

C{u) = / L{x,u^'^\x))dx, (11) 



Jn 

where Q, is an open, connected set in with smooth boundary 9f2 and L is a smooth function on 
the n-jet space (f2x]R^)^"^ A variational symmetry group of pT|) is a local group of transformations 
G acting o\i M d Q. xW with the following property: let fii be any subdomain with Vli C il, 
/ a smooth function on fii whose graph lies in Af , and suppose that for g G G, f :— g ■ f is a 
well-defined function on some rt[ G fl, then 

/ L{x\w'-''^f'{x'))dx' = ( i(x,pr(")/(a;))dx. (12) 
Jn{ Jni 

Heuristically, therefore, a variational symmetry group is a transformation group that does not 
change the variational integral. Since in the generalized setting we will exclusively deal with 
projectable group actions, let us investigate in the case where G is of the form [x, u) i— *■ 
(S^(x), ^,,(x, m)), S, smooth, a (w.l.o.g. orientation preserving) diffeomorphism for each rj. 
Then x' — S^(x) is a coordinate transformation and we can rewrite (|12p as 

i(S„(x),pr(")/'(S^(:E)))detDS,(a;)da;= / i(x, pr(")/(a;)) da: (13) 
1 Jni 

(with f — x' ^ 5'^(S_^(x'), /'(S_^(a;')))), for each fii as above and rj sufficiently small. 
Our strategy basically will be to turn (|13p into a definition in the generalized setting. Thus we want 
to allow for G a generalized projectable group action $ e tJ^(M x Rp+9)p+9 of the form (77, (x, u)) 
I— !■ (S^(x), \I',,(x, m)). Since contrary to the smooth setting, we cannot argue by continuity in r/ 
(in fact, we explicitly want to allow e.g. jumps), in order to make (jl3p well-defined for general 
L e Grii^ X R?)(")), we henceforth assume that ^ G ^[(— f7oi??o) x i^,^^] for some 770 > 0. This 
guarantees that for each f2i with Hi C there exists some 771 G (0, 770) such that 2^(77, a;) remains 
in a fixed compact subset of SI for all 77 < 771 and e sufficiently small. Therefore, the following 
definition makes sense: 

Definition 4.1 Let $ = (77,(0;,?/)) 1-^ (S^(a;), 5',,(a;, u)) G Qr{^ x W+<i)p+i he a generalized 
(projectable) group action with '^\{-r)o,rio)>^n G 5[(~%i%) x S7, f2] for some r]Q > 0. $ is called a 
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variational symmetry group of the functional with C £ ^/^-((rixIR')'^")) if for every Qi (- fli G Q 
and every f £ Q{fli)'^ we have (with f = x' ^ (a;'), : 

L(S,(x),pr(")/'(S^(x)))detDS^(a;)da;= / L{x,^y''^^ f{x)) dx mg{{-rn,m)) (14) 

whenever r]i > is such that remains in a fixed compact subset of fl for rj < rji and e 

small. 

Remark 4.2 (Inheritance properties of variational symmetries for smooth functionals) 
Consider a smooth variational problem (([3]) with L smooth) and a smooth projectable variational 
symmetry $. Then since composition of smooth functions and c-boundcd generalized functions is 
well-defined, holds for any / G C/[rii,R'^], since it holds on the level of representatives. The 
analogous statement is true for / g Q{0,i)'^ if we assume that both L and G are slowly increasing, 
uniformly for x in compact sets (to make the respective compositions well-defined). 



Remark 4.3 As was mentioned in Remark I3.2[ an alternative admissible set over which the 
variational problem C can be considered is the set of all c-bounded generalized functions from ft 
to W such that all partial derivatives up to order n are also c-bounded. In this case, in order for the 
transformed function /' to be well-defined we have to suppose that e ^[R x Rp+«, (Mp+'?)(")] 
(as can readily be verified by direct calculation using the chain rule). 

Our next aim is to derive an infinitesimal criterion for variational symmetries. In its proof we will 
make use of the following lemma. 

Lemma 4.4 Let S € Gt{^ x be a Q-complete group of transformations acting on R x Rp 

with infinitesimal generator ^. Then 

-^(det J(ry, x)) = det J(ry, x) ■ Div^(S„(a;)), 
dr] 

where J(r],x) is the Jacobian matrix ofE. with entries 

d 



= /t(s;(x-)) e g.(R X wy, i < i,j < p. 



and Div^ = X^iLi ^i^i ^•s total divergence of ^. 
Proof. 

Since the proof of the corresponding result in the smooth case given in [321 Theorem 4.12] can 
not directly be recovered in the generalized setting we give a direct computational argument. We 
have 

d " d 

— (det J(?7,a:)) = ^ ^ sgn(cr) Ji^(i) (?7, x) • . . . • — Jfe^(fc) (77, x) • . . . • J„<^(„) (?7, a;). (15) 



Here, 



d 9 " ^^fc 

-^Jka(k){l.x) = ^^(e''(2^(a:))) = ^^(S^(x))Jz,(fc)(r/,a;) 



Inserting into (jl5p we obtain: 

— {detJ{T],x)) = X! X! sgn(cr) Ji^(i)(?7,a::) • . . . • Jfe^(fc)(?7,x) • . . . • J„^(„)(?7,x)--^(S,,(a;)) 



Q^k \ 



l^k 

det J{r],x) ■ Div^(S^(x)) + 0. 



□ 
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Theorem 4.5 (Infinitesimal criterion) Let $ € Gt(^ x RP+'?)f+'' be a generalized, projectable, 
strictly Q -complete group action of the form given in Definition \4-l\ Then $ is a variational 
symmetry group of the functional (0) if and only if 

pr("V(L)+i-Dive = in X M«)(")), (16) 

for the infinitesimal generator 

P g Q g 

2—1 Q — 1 

Proof. Suppose first that $ is a variational symmetry of and let / G Q{fli)'', where fti C 

Cli C il. Then there exists some 771 > such that for all fj E (— ?7i, 771)^ and all C fii, by (I14p 
we have 

/ L(S^(a;),pr(")/'(Si;(x)))detDS^(a;)da; = / L(.t, pr(")/(a;)) int. (17) 

Since this holds for all open subsets il2 of ili, we may apply [44l Theorem 1] to conclude that for 
each rj as above the integrands in (fT7|) agree in Q{{Qi x M'')'")), i.e., 

L(S^(a;),pr(")/'(S,^(x)))det(J^(:r)) = L(:r,pr(")/(a:)) in g((r!i x M«)(")). (18) 

Suppose now that (ijM^"^) is any compactly supported gener alized point in (fli x M«)("). Then 
we may choose for the components of / G Q{Qi)'' suitable Taylor polynomials with generalized 
coefficients such that (S.u'"^) = pr(")/(2;). Therefore by definition of the prolonged group action 
we have 

L(pr(")$^(£, ?i("))) det(J^(£)) = i(£, «<")) . 

Since this holds for any given (a;,{t*^"-') £ ((ili x M*)("))~, by the point value characterization of 
generalized functions we conclude that 

L(pr(")$^(x, «("))) det(J,^(a;)) =L(a;,u(")) in ^^((-r/i, 771) x {fli x «<?)(")) (19) 

Differentiating this with respect to 77 yields (by Lemma [4.41) : 

(pr(")v(i) + iDivO • dct(J„(.T)) = in ^^((-7/1, 7^1) x (Q^ x R*)(")). (20) 

If we set ?7 = 0, we obtain (|16p on fli. Since fli as above was arbitrary, the claim follows on all of 

n. 

Conversely, if p6|) holds in G{{fl x R'')^")) for the infinitesimal generator ^ of $ then (j20p is 
automatically satisfied for each pair (17i,77i) as above. Thus we may integrate ([20]) from to 
?7 to obtain Evaluating at the graph of the generalized function / e Q{i}i)'^ (cf. (THl 

Definition 4.3.9]), again by the point value characterization in Q we obtain (jlSp . and, therefore, 
(|17p . from which we conclude that $ is a generalized variational symmetry group of £. □ 

As a preparation for Nother's theorem, we briefly recall the prolongation formula in Q from [161 
Theorem 4.3.17] (based on gl Theorem 2.36]): if v = T,LiC{x)-£t + ELi V'a(a;, u) gf^ is a 
tj — 71— complete generalized vector field with corresponding projectable group action $ on R^ x R"^ 
then 

pr(")v = v + ^^V'o(a;,"'"^)9„5, J ^ [ji, ■ ■ ■ ,Jk),l < Jk < P,l < k < n, 

a=l J 

where 

^i(x,zi(")) = Dj{^^ -Y.e<) + jZ^'u% e GriiW X R«)(")). 

i=l i=l 
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As in '55', (2.48)], we may introduce the characteristics of v as 

p 

g„(x, Va(x, u)-J2 Cut G Sr((Mf X 1 < a < g, 

i=l 

and Q[x, (Qi, . . . , Q,){x, m^)). Then = DjQ^ + Yll=i fwj,, and 

C(=l J J i=l a=l J ■■' 

= pr("VQ+f]fA, 

i=l 

with Di the i-th total derivative, and 

a— 1 a— 1 J 

Nother's theorem provides a relation between variational symmetry groups and conservation laws 
for the corresponding system of Euler-Lagrange equations. Recall that a conservation law for a 
system of differential equations A(x, u^")) = is a divergence expression Div P — Q which vanishes 
for all solutions u — f{x) of A (P is a p-tuple of generalized functions Pi{x, u'"') £ tJ((R^'+* 
l<i<p). 

Theorem 4.6 (Mother's theorem) Let ^ be a projectable, strictly Q-complete variational symmetry 
group of the variational problem and let 

i—l a— 1 

he the infinitesimal generator of $ with characteristics Qa o,s above. Then Q ■ E{L) is a conser- 
vation law for the Euler-Lagrange equations E{L) = 0, i.e., there exists a p -tuple = 
(Pi,...,Pp)(a:,M(")) such that 

DivP = Q- E{L) = ^Q^E^{L) in gr{{^ x 

Proof. The proof of the classical case ([42l Theorem 4.29]) carries over verbatim to our setting. □ 

For the practically most relevant case of first order variational problems, the explicit form of P is 
(cf. ga Cor. 4.30]): 

OL — 1 '' Q — 1 j — 1 * 

5 Applications 

5.1 Geodesies as solutions of a variational problem 

Let M be a smooth g-dimensional manifold equipped with a generalized Riemannian metric 
g e G2{M) (cf. [Ml [IS]). We are looking for curves c e t/[[0,T],M] which minimize the length 
functional ^ ^ 

L{c):^ [ \c{t)\dt= [ ^g{c{t),c{t))dt 
Jo Jo 
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As in the smooth setting (of., e.g., [H]), since L(c) < \J2T^ E(c) it suffices to find niinimizers of 
the energy functional 

E{c) := i £ |c(t)p dt^^ j\{c{t),c{t)) dt. (21) 

In terms of a local coordinate system {x^ , . . . , x"), the Euler-Lagrange equations for E read 

d'^c': y-^ .. del del , , , 

with rf^ the ChristofFel symbols of 5, 

gy = 5(g|T,g|7), and (g*-') = (5^^)^^. These are the usual geodesic equations for the metric 
5, whose derivation from (|2ip in the Colombeau setting, based on Lemma 13.31 is completely 
analogous to that in the smooth case (^21j). ((22|l is a system of ordinary differential equations for 
the components of c in tj([0, T]). 

Analogously, for pseudo-Riemannian generalized metrics, (f22|) is derived as the Euler-Lagrange 
equation for E by fixing the causal character (cf. |32j ) of c. 

As a concrete example let us consider the case of impulsive pp-waves (plane fronted gravitational 
waves with parallel rays) ([2l 1451 [27] ). Here, the metric is given in the form 

g = /(x, y)S{u)du^ — dudv + dx^ + dy^ 

with / smooth and 6 the Dirac measure. Using a generic regularization of 5 to embed g into the 
Colombeau algebra, one can show that the corresponding system of geodesic equations is uniquely 
solvable in C/[R]. Moreover, this unique solution is associated to the physically expected solution 
(piecewise straight geodesies which are broken and refracted at the shock hypersurface u = 0). 
We refer to the above papers or also to [TSl Sec. 5.3] for a detailed exposition. 



5.2 Examples from mechanics 

Denote by x{t) the position of a particle of mass 1 in M" that moves in a potential V{x). Accord- 
ing to Hamilton's least action principle, the particle follows the path that minimizes the action 
functional ^ 

£{x) = / L{x{t),x{t))dt (23) 



among all trajectories y(t) such that y{a) — x{a), y{b) — x{b) for whatever points of time a < b. 
Here the Lagrangian is of the form 

L{t,x,x)^^\if-V{x) (24) 

and given by the difference of the kinetic and the potential energy. It is clear that time translations 
form a variational symmetry for the functional ((23l) . It follows from Nother's theorem that the 
total energy 

P^l\i\^ + V{x) 

is conserved. Particle systems are of interest here, because our approach offers the possibility 
to allow singular potentials, for example, delta function potentials. Singular potentials may be 
modelled by elements V G tJr(R")- If a; e (^(M))" minimizes the action functional on every 
interval [a, 6], then it follows from Theorem 13.41 that x satisfies the Euler-Lagrange-equations 

x{t) +VV{x) ^0. (25) 
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It has been shown in [16] that (j25p . completed by initial data x{tQ) = xq, i(to) = yo bas a 
solution X e (^(M))", provided V is either of i°°-type or nonnegative. In addition, total energy 
is conserved along the generalized trajectories. 

Uniqueness of the solution to in (tJ(M))" can be obtained when W^V is of L°°-log-type. As 
in the classical case, the second variation of (|23p is rarely positive so that uniqueness cannot be 
inferred using the variational principle alone, in general. 

We now turn to a more explicit example - a classical particle in a delta function potential. For- 
mally, it is described by the Lagrangian L{t I — S(x). The delta function potential 
is considered as acting as a barrier so that the particle trajectories exhibit pure reflection at the 
origin. To make things more precise, let us call a generalized function D G ^?r(K) a strict delta 
function if it has a representative {pe)e G fT(R) such that 

(i) supppe {0} for e 0. 

(u) lime^o / Pe{x)dx = 1. 

(iii) J |/9£(a;)| dx is bounded uniformly in e. 

D is called a model delta function if it has a representative ((/?e)e of the form '^e(x) = ), 
where G / '^{x) dx = 1. 

Example 5.1 (Particle in a delta function potential) The Euler-Lagrange equations are 

x{t) + D'{x{t)) = 

x(0) ^ xo (26) 
i(0) = ya 

where D E Gt(^) is a strict or model delta function and Xo,yo € R. We have the following result: 

(i) Problem ^ has a solution xeg{R). 

(ii) If D" is of i°°-log-type, the solution is unique in ^/(M). 

(iii) If D is a model delta function and xq = XQ^yo = j/o are elements of M with xq ^ 0, then the 
solution X e t/(M) of ([26| is associated with the function 1 1-^ sign(xo)|a;o + tyol- 

The pure reflection picture does not necessarily apply when strict delta functions are employed. 
In general, it is possible that certain trajectories terminate at the origin, so that the particle may 
be trapped there after finite time. More details and proofs of these results can be found in [16] . 

Example 5.2 (Planar motion in a central force field) Introducing polar coordinates x — rcosip, 
y ~ r sin ip in the plane, the motion of a particle of mass m in the time interval [a, b] is a minimizer 
of the action functional ^ 

£(r,(^) = / L{r{t),r{t),^{t),^{t))dt 



with 

L(r, f, (f, Lp) = -m{f'^ + r^ip^) - V{r). (27) 

Here the potential energy V{r) is assumed to be an element of C/t-(R). Such a situation could 
arise from distributional potentials as in the previous example, or from regularizing classical, non- 
smooth potentials like the gravitational potential — 1/r. It is clear that the infinitesimal generator 
of the rotation 

- A 

dip 

determines a variational symmetry. According to Nother's theorem fTheorem l4.6l) the correspond- 
ing conserved quantity is angular momentum 

F = — — — mr If, 
op> 

and this holds for generalized solutions r,Lp ^ t/[a, b] as well. 
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5.3 Examples from elastostatics 



We begin with an example from classical elastostatics. To keep matters simple, we do not treat the 
full three-dimensional problem but rather focus on a model problem that describes, e. g., vertical 
displacements of membranes (n = 2) or strings {n = 1). 

Example 5.3 (Standard linear elasticity) Let fl C M" be a bounded, open and connected set with 
smooth boundary dil, the reference configuration of an elastic body. We denote by u : R 
the displacement. For simplicity, we address the Dirichlet problem = uq only. The static 
solution minimizes the functional 

C{u)^l- I {a{x)\Wu{x)\^ + f3{x)u{x)^) dx ~ j j{x)u{x)dx 
2 Jn Jn 

under all displacements that satisfy the boundary condition. Here a{x) depends on the elastic 
properties of the material (modulus of elasticity and Poisson's ratio), f3{x) is a potential that 
could arise by an elastic constraint, for example, and j{x) is the body force. In our setting, we 
may take a, (3 and 7 as generalized functions. In this way, singular potentials and degenerate 
coefficients a{x) can be modelled. A delta potential P{x) — S{x — xq) would arise if a spring is 
attached to a membrane at point xq € f2. If the coefficient a{x) vanishes in a subregion of f2, 
the classical problem looses wellposedness, but the generalized problem might still be wellposed if 
a{x) is modelled by a generalized function associated with zero, but still invertible. 
As in Example 13.81 we make the following assumptions: a,/3,7 € Q{^), a > and invertible on 
f3 > 0, uo (z G{dQ). The admissible set and the admissible variations, respectively, are given by 

A={ue g{n) : u\on = uo}, Aq ^ {u E g{n) ■. u\on 0}. 

Observe that C is the quadratic functional 

^(u) = ^aiu,u) - f{u) 
corresponding to the the bilinear form 

a{u, v) = [a{x)\i' u{x) ■ \7v{x) + (3{x)u{x)v{x)) dx 
Jn 

and the linear form 



The Euler-Lagrange equation is 

- div [a{x)Vu{x)) + P{x)u{x) = 7(0;), u\gn = uq. (28) 

The first and second variation was calculated in Example 13.81 and it was shown that C(u) admits 
a unique minimizer. It is clear that u € A solves the Euler-Lagrange equation (|28p in g{Q) if and 
only if it solves the variational equation ([8]). In this way, we have proved that solutions to the 
Euler-Lagrange equation are unique. 

In case the coefficients a, (3, 7 are classical smooth functions, a proof of existence can be obtained 
as follows: Fixing a representative uq^ of uq, equation (j28p has a unique classical smooth solution 
Me, e G (0,1], cf. [H], Chap. 6]. The closed graph theorem yields continuous dependence on 
the boundary data in the C°°-seminorms. This shows that {ue)e is moderate and thus defines a 
solution u S g{i}) to (|28|) . An existence result with generalized coefficients a € G{^), /3 = 7 = 
has been obtained in [43] . Alternatively, the solution can be constructed by direct methods of the 
calculus of variations, as have been worked out in the Colombeau setting in jllj . 



/(") 



"f{x)u(x) dx. 
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Example 5.4 (Euler-Bernoulli beam with discontinuous coefficient) The equilibrium equation 
for the deflection w{x), x € of an Euler-Bernoulli beam under constant axial force /3 and 

distributed load j{x) is 

cP / , ^(fw(x)\ ^cPw(x) , . „ r„ „i 

Here £ is the length of the beam. The coefficient a is given by the flexural stiffness EI [E the 
modulus of elasticity, / the area moment of inertia) and may depend on the position x along the 
beam, in general. Appropriate boundary conditions for a beam freely supported at both ends are 

w{Q) = 0, w{e) = 0, -^wiO) = 0, -^wit) = 0. (30) 

We are particularly interested in the following cases: 

(a) the coefficient a has jump discontinuities (non-smooth case); 

(b) the coefficient a vanishes in certain points (degenerate case); 

(c) the coefficient a is given by a path of positive noise. 

Case (a) is a simple model of abrupt change of material properties (see e.g. [20l |47]); case (b) 
arises for example, if the beam contains a joint in its interior. Classically, these problems are 
split into adjacent subproblems with transition conditions. A global formulation in the classical 
setting - particularly in case (b) - is impossible, because the classical ellipticity condition would be 
violated. Case (c) is of relevance in stochastic structural mechanics. Since a has to be nonnegative, 
a lognormally distributed noise process would be a suitable model for random disturbances. Such 
a process does not exist in the classical theory of stochastic processes, but it does exist in the 
sense of a Wick exponential [H] or as a Colombeau process [41] . In all three cases the Colombeau 
setting lends itself as a framework that provides a rigorous global solution concept. 
We observe that ((29|) is the Euler-Lagrange equation corresponding to the quadratic functional 

1 

C{u)^- {a{x)\u"{x)\'^ - (3\u'{x)\^)dx- -/{x)u{x)dx. 
2 Jo Jo 

The second variation is given by 

1 

S^Ciu,v)^- {a{x)\v"{x)\^ - I3\v'{x)\'')dx (31) 
2 Jo 

We assume that there is an invertible element ao > in M such that a > ag in C/[0,^]. Clearly, 
the second variation ((3T|) is positive definite if the axial force f3 is negative. We are going to show 
that pTjl is positive definite even for small, positive axial forces /3, below the first eigenvalue of 
equation To see this, let w e C°°[0,£], w{0) = w{£) = 0. Then 

w{x)=--J (i-y)w"(y)dy+ I {x ~ y)w" {y) dy . 
From there we derive the inequality 







II"'IIl2(Q.^) < C'lh"llL2(0/) 

for some constant C > and all smooth w as above. Further, 



w' {x)\'^ dx — I w" [x)w{x) dx 
Jo 



- -\[w'i^)\''dx-\j\w{x)\'dx 
> -1(1 + C) j\w"{x)\'dx. 
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From this inequality it follows that (|3ip is positive definite if ap — ^(1 + C)/3 > in Q[0, £\. If this 
expression is invertible, then C admits a unique minimizer (Prop. [321) i which is obtained as the 
solution to the Euler-Lagrange equation (|29p . 

As a concrete application, we are going to work out a global formulation of an Euler-BernouUi 
beam with a joint, say at half-length x = £/2. At a joint, the flexural stiffness , i.e. the coefficient 
a, vanishes. In the Colombeau setting, this can be modelled as follows. Take a smooth unimodal 
function vanishing for \x\ > 1 and equal to 1 for |a;| < 1/2, < ip{x) < 1 and let ~ define a 
nonnegative invertible element of R. Put 



aAx)^a[l-{l-h,)^{^^^)) 



where a e M is some positive constant. Clearly, {ae)e determines a nonnegative invertible element 
of CJ([0,^] which is infinitesimally small near x = £/2 and otherwise equal to a. We consider the 
Euler-BernouUi beam without axial force 



dx^ 



(ae(x)^(x)) -7(x)=0, xe[0,£] (32) 



where 7 G C°°[0,£] denotes the distributed load. A representative Ug of the solution in tj[0,£] is 
readily calculated as 



Jo aeiz) £ Jo Jo ae{z) 

where M{x), minus the bending moment, is the solution to M"{x) — 7(2;), M(0) = M{£?) = 0. 
Now, 

^M(z)^ \ r.^dz, 0<y<^/2-e, 
dz — 



a.{z) ^ M£) + + /,%^, ^ dz, £/2 + e<y<£ 



with 



Je/2-e Cte{z) 

If we choose in such a way that 

dz edy 



D (33) 



/2^e l-(^-h,)^P{^) 7_i 1 ~ (1 - /le)V'(y) 
in R as e ^ 0, we will have that 

fyUiz)^ \ lo'^dz, 0<y<£/2, 
lim / — dz = < 

Jo ^e{z) \ fj^ M£) + «i5M(|) + M£) d,, £l2<y<£ 

and so Uf,(x) will converge to a continuous, piecewise smooth function u(x). This limit - the asso- 
ciated distribution - describes the displacement curve of the beam. Larger values of D correspond 
to larger loss of stiffness at a; = ^/2 and thus to larger displacements of center of the beam. Failure 
of the beam can be modelled by letting Z? = 00 in p3)) . 

We remark that the case of a discontinuous fiexural stiffness and nonzero axial force has been 
treated in [SOj. Finally, in the stochastic case, positive noise on [0,^] can be defined in the 
Colombeau sense as a lognormally distributed process with mean value 1 and variance given by 
exp(||pe||^2(]{)) — 1, where is a strict delta function as in Subsection 15.21 The paths of this 
generalized random process are nonnegative elements of the Colombeau algebra ^/[O,^] and thus 
may serve as describing highly random behavior of the flexural stiffness a. We refer to [41] for 
details on and explanations of the positive noise model. 
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Example 5.5 (Rods with generalized stress-strain relation; hard rods) Let u{x) be the displace- 
ment of a rod of length £ of cross-sectional area one, subject to the body force (density) f{x), 
< X < £. The balance law is 

(T'{x)+f{x)^0 

where cr denotes stress. Letting e — ^ the strain, assume a stress-strain relation (constitutive 
law) of the form a = (?(e). If the rod is clamped at one end and free at the other, the displacement 
is the solution to the problem 



If the constitutive law has a potential, g{y) — —G'{y), equation p4|) is the Euler-Lagrange equation 
of the functional ^ 

C{u) = (g{u'{x)) + f{x)u{x)^dx. (35) 

We allow generalized potentials G G tJr[0,^] and forces / e ^/[0,£]. The admissible set and 
admissible variations are 

A = Aq = {ucg[Q,l] :u(0) = 0}. 
Assume that u ^ A minimizes the functional (j35p . The first variation is 

&C{u-v) = j g{u'{x))v'{x) + f{x)v{x)^dx 

[g{u'{x))'v{x) + f{x)v{x)) dx + g{u'{l))v{e). 



If SC{u; v) — for all v £ Aq, we have in particular that 

g{u'{x)y + f{x)=0 

in Q{0,£), by using Lemma 13.31 In order to show that g{u'{£)) = we need an adaptation of 
the proof of this lemma. Thus assume that g{u'{£)) 7^ in M. We can find {re)s € Ma/ and a 
subsequence ^ such that ge^ (u^^ {£)) r^^ = 1 for all fc e N. On the other hand, there is iV > 
such that 

sup \geKix)Y + fe{x)\ = 0(6"^) as £ 0. 
o<x<e 

Similar to the argument in Lemma 13.31 we can find an element v £ G[0,£] with support in {0,£] 
such that {£) = r^^ and 



{g^(u'^{x)) + f^{x))ve{x)dx 



1 

< - 
- 2 



for all e, contradicting the hypothesis that S£{u; v) = for all v £ Aa- We arrive at the conclusion 
that a minimizer of (|35p in ^ is a solution to the Euler-Lagrange equation 

In linear elasticity, the stress-strain relation is <t = Ee with the modulus of elasticity E. The 
larger E, the harder is the rod. In the limiting case E' — > oo, the rod becomes inextensible. We 
can model such a rod in the Colombeau framework by introducing the potential G £ Gt{^) by 
Ge{y) = — \y- The Euler-Lagrange equation becomes 

-u'^[x) -t- A(x) = 0, u,(0) = 0, u'M = 

£ 

and has the solution ^ 

Us{x) =e / fe{z)dzdy. 



For bounded body force, the solution is associated with zero, thus no extension of the rod takes 
place, as anticipated. We refer to |33j for the corresponding model in convex analysis. 
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5.4 Nonlinear wave equations with singular potential 



Example 5.6 (Wave equation with delta potential) In this subsection, we discuss the one-dimen- 
sional wave equation 

at2 " 9^""^^' 9^^^"^' " ° ^^^^ 

with potential V{x^u). The higher dimensional case can be treated similarly. Classical examples 
of the potential are V{u) = + ^u'^, m, fc > 0, leading to the cubic Klein-Gordon equation, 
or V{u) — — cosw, leading to the Sine-Gordon equation. However, much more singular potentials 
are in use. For example, the equation 



-Q^u{x, t) - -^u{x, t) -I- 5{x - xo)F{u{x, t j) = 



(37) 



describes the vibrations of a string with a nonlinear spring with restoring force F(u) attached at 
the point x = xq, see [22j . An even more singular potential. 



u{x, t) 



u{x, t) — 6{u{x, t) — Uk) — 



fc=i 



where {ui, . . . ,u„} is a finite subset of K, has been studied in [5]. This equation arises from a 
piecewise constant approximation of a smooth potential - in the quoted reference, the Ginzburg- 
Landau potential. Note that in this example, the sought after solution has to be inserted in a sum 
of Dirac measures. All these cases can be subsumed by letting the potential V{x, u) belong to the 
Colombeau algebra Gt{^ x R). The Lagrangian functional corresponding to ([36|) is 



Ciu) = 



d 2 
^u{x^t) 



d 2 



V{x, u{x, t)) j dxdt 



The element u G t/(R^) that makes C{u) stationary on every rectangle [a,b] x [c,d] produces a 
solution of the Euler-Lagrange equation ([551) ■ 

Existence and uniqueness of solutions in Q{W^) to the Cauchy problem for equation p6p have been 
demonstrated under various conditions on the potential V, see e.g. j37[ 134] . Solutions to equation 
([57)1 in the Colombeau algebra t/(M^) and their associated distributions have been computed in 
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